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It has been recently found that the equations of motion of several semiclassical systems must take 
into account terms arising from Berry phases contributions. Those terms are responsible for the 
spin Hall effect in semiconductor as well as the Magnus effect of light propagating in inhomogeneous 
media. Intensive ongoing research on this subject seems to indicate that a broad class of quantum 
systems may be affected by Berry phase terms. It is therefore important to find a general procedure 
allowing for the determination of semiclassical Hamiltonian with Berry Phase corrections. This 
article presents a general diagonalization method at order h, for a large class of quantum Hamiltonians 
directly inducing Berry phase corrections. As a consequence, Berry phase terms on both coordinates 
and momentum operators naturally arise during the diagonalization procedure. This leads to new 
equations of motion for a wide class of semiclassical system. As physical applications we consider 
C ' here a Dirac particle in an electromagnetic or static gravitational field, and the propagation of a 

^ ' Bloch electrons in an external electromagnetic field. 
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^ ■ I. INTRODUCTION 

> 

' Since the seminal work of Berry 1] , the notion of Berry phase has found several applications in branches of quantum 
T-H : physics such as atomic and molecular physics, optics and gauge theories. Most studies consider the geometric phase 
CO that a wave function acquires when a quantum mechanical system has an adiabatic evolution. Yet, the Berry phase 
in momentum space has recently found unexpected applications in the topic of spintronics. Such a term may indeed 
be responsible for a transverse dissipationless spin-current in semiconductors in the presence of electric fields [5]. 
This effect is a particular case of the Spin-Hall effect which is now predicted and observed in many different physical 
situations and can be interpreted at the semiclassical level as due to the influence of Berry connections on semiclassical 
equations of motion of spinning particles, like electrons in electric fsl or magnetic fields In the above cited examples, 
Oh, the semiclassical equations of motion where derived from an approximate semiclassical Hamiltonian in a representation 
^ where this latter is diagonal, ft was then shown that a nonconimutative geometry, originating from the presence of a 
Berry phase which turns out to be a spin-orbit coupling, underlies the semiclassical dynamics. Spin-orbit contributions 
on the propagation of light have also been the focus of several other works [1, @| and have led to a generalization 
of geometric optics called geometric spinoptics [7] . 
^ Semiclassical methods play a very important role in solid state physics too, in studying the dynamics of electrons 

\ to account for the various properties of metals, semiconductors and insulators [1]. In a series of papers Q (see also 
[lo|), a new set of semiclassical equations with a Berry phase correction was proposed to account for the semiclassical 
dynamics of electrons in magnetic Bloch bands (in the usual one band approximation). These equations were derived 
by considering a wave packet in a band and using a time-dependent variational principle in a Lagrangian formulation. 
The derivation of a semiclassical Hamiltonian was shown to lead to difficulties in the presence of Berry phase terms 
0. The apparent non-canonical character of the equations of motion with Berry phase corrections led the authors 
of [ni to conclude that the naive phase space volume is not conserved in the presence of a Berry phase and a 
magnetic field. There is nevertheless an invariant measure as found by themselves, so the Liouville theorem is not 
really violated in the end. This invariant measure is actually a general result of the well-established theory of non- 
canonical Hamiltonian dynamics, as pointed out by a number of authors [T^ . [Tsl . [l^ . However it is only in fl^ 
that the non-canonical Hamiltonian formulation of a semiclassical electron in magnetic Bloch bands has been fully 
derived. This Hamiltonian approach allows deriving rigorously the semi classical equations of motion, including 
explicitly the role of the Berry curvature and showed many similarities with the description of a Dirac electron in 
an electromagnetic field as in Q . The common feature of the Hamiltonian formulations discussed above is that a 
noncommutative geometry underlies the algebraic structure of both coordinates and momenta. Actually, a Berry 
phase contribution to the dynamical operators stems from the representation where the kinetic energy is diagonal (for 
instance Foldy-Wouthuysen or Bloch representation). In this representation the physical coordinate and momentum 
become noncommutative operators. 
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The previous discussion shows that Berry phase terms could be present in semiclassical equations of motion of several 
physical systems ranging from electrons in vacuum, in solid or in semiconductor to photons in inhomogeneous media, 
with potential application in the field of spintronics and spinoptics. This in turn called for a general semiclassical 
Hamiltonian formalism from which semiclassical dynamics of a quantum system can be derived. This paper presents a 
general method of diagonalization at order h for a quantum mechanical Hamiltonian presenting bands structure, like 
for electrons in a periodic potential or for a Dirac (massive or massless) like-Hamiltonian. Starting with a Hamiltonian 
depending only on an invariant momentum P and whose diagonalization is known, we can introduce a dependence 
in the variable R and diagonalize the Hamiltonian in four steps (discussed in the text). During this process of 
diagonalization, we show that both position and momentum operators acquire a Berry-phase contribution making 
both the coordinate and momentum algebra noncommutative. As physical applications and to check to validity of our 
method we further consider the case of Dirac particle in an electromagnetic field and compare with the semiclassical 
diagonalization given in Q . We also consider the case of a Dirac particle in a symmetric static gravitational field (the 
asymmetric case is studied in ref. [l6|) and compare with the articles [13, Ell- As a last apphcation, the reader will 
find the details of the diagonalization sketched in [l^.for the propagation of a Bloch electron (spinless) in an external 
electromagnetic field. These various physical applications show that our semiclassical Hamiltonian diagonalization 
approach is potentially promising since it unifies several apparently unrelated problematics in one formalism. 

The paper is organized as follows. In section II we develop our formalism in the case of a general Hamiltonian 
which has an energy bands structure. We then derive the very general equations of motion in this case. Section III is 
devoted to the application of our method to the case of the Dirac Hamiltonian in an electromagnetic field but in a fiat 
space, and to the diagonalization of the Dirac Hamiltonian in a symmetric static gravitational field, allowing us to 
check the validity of our method. Section IV retrieves the equations of motion for an electron in a periodic potential 
within our general set up. Section V is for the conclusion. 

II. A GENERAL PROCESS OF SEMICLASSICAL DIAGONALIZATION. 

In this section we present a method to diagonalize at the semiclassical order (h) a quantum mechanical system 
whose state space is a tensor product {T^^) ® V with V some internal space. In other words, the Hamiltonian of 
this system can be written as a matrix Hq (P, R) of size dim V whose elements are operators depending on a couple 
of canonical variables P and R. The archetype example is the Dirac Hamiltonian with V = C**, but in appendix 2 we 
show how a spinless electron in a periodic potential fits also in this set up. By diagonalization, we always mean here 
a unitary transformation setting the Hamiltonian in a diagonal matrix form, the diagonal elements being operators 
depending on P and R. That is, we do not aim at finding the eigenvalues, but rather to derive the Band Hamiltonians, 
that are usually relevant for the semiclassical dynamics. 

More precisely, our strategy will be to solve the exact diagonalization for Hq (P, R) at order h, when the diagonal- 
ization of a fictitious Hamiltonian Hq (P,r) is known for a parameter r (replacing R) which is supposed to commute 
with P (for instance this may be the Foldy-Wouthuysen transformation [19| for a free Dirac particle). We show how 
to compute the quantum corrections (at leading order in h) that were neglected during this formal diagonalization 
(where position and momenta where considered as commuting quantities). The idea behind this procedure is that it 
is much easier to solve the diagonalization for Hq (P,r), as seen in our applications, and only then turn to Hq (P, R). 

A. Preliminary : Products of operators series. 

To develop our process of diagonalization, the semiclassical expression of products of symmetrized expressions 
S (P, R) depending on the canonical couple of variables P and R is required. These expressions are also assumed to 
have series expansions in P and R whose coefficients can be of a matrix form (this last assumption allowing to deal 
with Dirac Hamiltonians). Let us consider two such expressions 5*1 (P, R) and ^2 (P, R), supposed to be symmetrized 
in P and R. By symmetrization, we mean that each expression has been written in a form where all the powers 
of P have been put half on the left and half on the right of the expression. Our aim is now to write the product 

(P,R) S2 (P,R) as a symmetric expression in terms of P and R. This is easy to realize at order h, since in that 
case, pushing half of the powers of P in ^2 (P, R) on the left and half the powers of 5*1 (P, R) on the right is equivalent 
to the computation of some commutators. One can easily see that at order h 

Si (P, R) ^2 (P, R) - SymiSi (P, R) ^2 (P, R)) + ^MsymV^, Vp, [Si (P, R) ^2 (P, R)] (1) 
where Sy'm{Si (P, R) 6*2 (P, R)) is the symmetrized expression for the product, and Asym is defined through 



Asym^n^Vpi [Si (P, R) ^2 (P, R)] = V^^^i (P, R) VpiS2 (P, R) - Vp,^i (P, R) Vp,52 (P, R) (2) 
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This formula can be easily generalized to an arbitrary product of n terms, but for the sequel of this paper, only three 
will be needed 



where 



5i(P,R)52(P,R)53(P,R) = 52/m(5i(P,R)52(P,R)^3(P,R)) 

i 

+ 2 

i 

+ 2 



hAsymVR,Vp^ [Si (P, R) ^2 (P, R)] ^3 (P, R) 
hSi (P, R) AsymVfl, Vp, [^2 (P, R) ^3 (P, R)] 



--hAsymVR,Vp, [Si (P, R)s,(p,r) ^3 (P, R) 



AsyrnVR^pi [Si (P,R)s,(p,r) ^3 (P,R) 
[Vfl, ^1 (P, R)] S2 (P, R) Vpi Ss (P, R) - [Vpi Si (P, R)] ^2 (P, R) V^, ^3 (P, R) 



(3) 



(4) 



Let us stress again that all this identities are valid only at order h and that considering higher order corrections would 
of course induce more corrections. 



B. Diagonalization with a parameter r 

Let us consider a general Hamiltonian Hq (P,R) which admits a series expansion in P and R written here for 
convenience in a symmetrical form. To perform the semi classical diagonalization of this operator, we first consider in 
this subsection a fictitious Hamiltonian Hq (P,r) where a parameter r commuting with P has replaced the operator 
R. We further assume that Hq (P,r) is known to be exactly diagonalized through a matrix U (P,r) = U. We assume 
that U and U (P,r) Hq (P,r) (P,r) can be expanded in series of P (with positive or negative powers) as it will be 
the case in our applications, but in fact this assumption which is practical for our proofs could probably be relaxed. 

As an typical example we can consider the following kind of Dirac Hamiltonian 

Hq ^a.(P- A(r)) +/3m (5) 

where A(r) mimics a formal magnetic potential. The potential A(r) being R independent, it only shifts the momentum 
for each value of r. As a consequence, the usual Foldy Wouthuysen [19] transformation expressed in terms of P — A 
instead of P, diagonalizes the Dirac Hamiltonian exactly. 

Going back to the general case, the diagonalization of Hq ( P,r) will be written as 

e (P,?) = U (P,F) Hq (P,?) U+ (P,7) (6) 

where e (P,r) is a diagonal matrix. For later use, let us notice that we have U (P,r) PU^ (Pi?) = P- 

To gain some hints from our initial diagonalization, recall we assumed that the product e(P,r) = 

U (P,r) Hq (P,r) (P,?) can be expanded in series of monomial terms of the form : (^e]^ (r).P^^^^^ ... ^ej^" (r).Pj^'" j 

(*) where the e]' (r) are some matrices commuting with P, and P^^'' is the ij-th component of P {ii = 1, 2, 3) at some 
power fci, (this power should be written ki^^i but we will avoid an excess of notation here). As usual, the sums over 
the ii are implicit. Rearranging the series in powers of Pi, P2, P3 and given that r is a parameter, we can write the 
energy in a symmetrical form : 

e(P,?) = i^ 

^ X 

where ^Oi^i 2 3 ^i^') ^ given momentum power and As.(j) is a combination of the e;'(r). The X labels the multi 
index {Xi, i = 1,2, 3). 

The important consequence here is that the matrix Ax(r) is diagonal. As an example, think of the Dirac Hamil- 
tonian diagonalization, which involves some products a.Pa.P, rearranged as 

i [aiPiajPj + ajPja^Pi] = ^ [aiUj + a^a^] PiPj (8) 

and i [aiaj + a^ai] is diagonal, by the usual rules for Dirac matrices. 

Let us conclude this subsection by noting that the symmetrizations we performed is unnecessary here, but will be 
of a practical interest when dealing with the exact diagonalization. 



Ax(?) n p/^' 



n 



p 



=1,2,3 



(7) 
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Introducing the R dependence. The transformation Ansatz. 



We are now going to reintroduce R into Ho in order to diagonalize exactly Hq (P, R) at the h order. In the 
foUowing symmetrization in R and P wiU be assumed in all expressions. To find the diagonalization transformation 
for Hq (P, R), we will use the following method. First notice that the Hamiltonian Hq (P, R) is "almost" diagonalized, 
that is diagonalized at zeroth order in h, through the transformation U (P, R) (which is also symmetrized given our 
convention) . 

However, through the symmetrization process, the matrix U (P, R) does not remain unitary. As a consequence, we 
will rather consider a matrix U (P, R) + XU (P, R), where AT is a contribution of order h that ensures the unitarity 
of the transformation at order h (the factor J7 is a normalization that simplifies the subsequent expressions). 

The matrix X can be explicitly computed. Actually, from the unitary conditions 



and 



or equivalently 



(C/(P,R) + AC/) (C/+(P,R) + C/+A+) 

(C/+(P, R) + U+X+) (C/(P, R) + XU) 

f/(P,R)C/+(P,R) + A + A+ = 
C/+(P,R)C/(P,R) + C/+(A + A+)C/ = 



1 



= 1 



(9) 
(10) 

(11) 



To solve this equation, let us first notice that [/(P, R)f/+(P, R) ^ 1 since C/(P,R) is not unitary. The crucial point 
here and in the sequel of this paper is the computation of a product of expressions in which the R dependence has 
been introduced. To do so let us use the initial relation for the parameter r 



C/(P,?)C/+(P,?) =1 
expanded in a symmetric series (as we did for the Hamiltonian) 



C/(P,7)C/+(P,?) = i^ 



X 



^x(?) n + n 



, 1=1,2,3 



,4=1,2,3 



= 1. 



(12) 



(13) 



Therefore the series expansion J2x ^x(r) (^ni=i 2 3 ^/^') reduces to one constant term : the identity matrix. 
Now, going back to J7(P,R), we use the symmetrization formula ([T]) 



f/(P,R)[/+(P,R) 



Sx(R) 




Bx(R) 



--hAsymVR,Vpi [C/(P, R);7+(P, R)] 



(14) 



Since our result about the Sx(r) applies also to i?x(R-) after replacing r by R, one has : 

?7(P,R)L/+(P,R) = l + ^hAsymWR,Vp, [f/(P, R)f/+(P, R)] 



(15) 



where we have defined the ("non projected", see below) Berry connections (we use the terminology of ref. [20|) as 

= i?iC/(P,R)VpC/+(P,R) 



Ap = -i;iC/(P,R)VRC/+(P,R). 
that are non diagonal Hermitian matrices of order h. Similarly, we have also: 

C/+(P,R)C/(P,R) = l + ^hAsymWB,Vp, [[/+(P, R);7(P, R)] 

= 1-U+{P,R)^[Ap.,Ar,]U{P,R) 



(16) 



(17) 
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Let us note that, as can be checked easily, the Berry phases are Hermitian. So is ^ [^pi,^/?,]. Therefore we can 
solve our problem with 

X+ - ^[Api,Ar,] (18) 

Let us make an important remark at this point. Our choice for X is obviously not unique. Actually, it has been 
chosen to ensure the unitarity of the transformation and to obtain a transformation that reduces to the initial one 
when P and R do commute. We could thus add to X an expression like 5XU where 5X is anti-Hermitian. It is easy 
to see that the operator U (P,R) + XU + 5XU is still unitary. However, we will soon see that this non unicity is 
irrelevant and that our choice is sufficient to perform the diagonalization at order %. 



D. The quasidiagonalization 

We will now consider the following quasi-diagonalization transformation 

[([/ (P, R) + XU) Hn (P, R) (C/+ (P, R) + U+X+)] 

To compute this last expression, decompose it at the first order in h as 

U (P, R) Ho ((P, R)) U+ (P, R) + XUHo (P, R) U+ (P, R) + C/ (P, R) Hq (P, R) U+X' 
~ U ((P, R)) Ho ((P, R)) U+ ((P, R)) + Xe ((P, R)) + s (P, R) X+ 



(19) 



(20) 



Let us first have a look to U (P, R) Hq (P, R) [/+ (P, R) and consider it, as before, as a series of products of operators 
in P and R.If this two variables where commuting, we would recover the expansion e (P, R) given in the subsection 
B. But now, since, R does not commute with P, one has rather 



[/(P,R)7?o(P,R)t/+(P,R) = 



Ax(R 



commutators 




Ax(R) 



(21) 



the commutators appearing while pushing the momentum powers on the left or on the right. 

We can compute the first two terms of the right hand side by the same trick as before. Actually, by construction, the 
coefficients of the series expansion of Ax(R) in the variable R, are the same as the coefficients (which are diagonal) 
of the expansion of j4x(r) in the parameter r. As a consequence. 



Ax(R) n + n ^""m ^x(R) 



i 1=1,2,3 



, 1=1,2,3 



(22) 



is the series expansion of e (P, R), the powers of P being rejected symmetrically to the left and to the right. 



As an example, consider again the case of the Dirac Hamiltonian with an electromagnetic field. The free 
"Benchmark" case is (P) = P^, and given our conventions, replacing P by P — A(R) leads us to define 
£2 (P - A(R)) = P2 - A(R).P - P.A(R) + A2(R), which is simply the usual operator (P - A(R))^ By the same 
way, we obtain as a series expansion 



s (P - A(R)) = e (P) - - ( A(R). — +— .A(R) 



1 



1 



(23) 



The last term in the right hand side of ((2T|) involves just half the commutators obtained in pushing the momentum 
operators to the left or the right. As explained in the subsection A Eq.®, they are simply given by 



[commutators] = -hAsymV b.^ p' [U (P, R) Hq (P, R) C/+ (P, R) 



(24) 
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As a consequence, we can write 

[/ (P, R) Ho (P, R) C/+ (P, R) = e (P, R) + ^MsymVfl, Vp, [U (P, R) (P, R) U+ (P, R)] (25) 
A lengthy but straightforward computation presented in appendix 1 leads to 

f/(P,R)ifo(P,R)f/+(P,R) = £(P,R) + i[^fl,Vfl,£(P,R) + VH,e(P,R)^ijJ 

+ i [^p,Vp,e (P,R) + Vp,e(P,R) Api] 

[£ (P, R) , Api] Ar, + ^[e (P, R) , Ar,] Ap, 
+^[^fl,,^p.]£(P,R), (26) 

To end up with the quasi-diagonalization, wc have to add the expression 

X£(P,R)+£(P,R)X+ (27) 
Given the expression obtained previously for X, we have thus 

Xe (P,R)+e(P,R)X+ 
= -^[^fi,,^pi]£(P,R)-^£(P,R)[^iip^P!] (28) 

We can thus ultimately write the diagonalization process as : 

[{U (P, R) + X) Ho (P, R) {U+ (P, R) + X+)] = £ (P, R) + 1 [AR,VRie{P, R) + Vr^s (P, R) Ar,] 

+1 [^pi Vp,£ (P, R) + Vp,£ (P, R) Api] 

[£ (P, R) , Api] Ar, + ^ [£ (P, R) , Ar,] Api 

+^[[^fl,,^p.],£(P,R)] (29) 

Let us conclude this section by noting that our transformation is not a diagonalization at order U, since it in- 
cludes non diagonal contributions of order h through the Berry connections Ar^ and Api , justifying the name quasi- 
diagonalization. However, the next paragraph will show that these non diagonal terms are only an artifact. Actually, 
projecting our transformed Hamiltonian on the diagonal will in fact yield the genuine diagonalization. 



E. The "Exact" semiclassical diagonalization 

As mentioned in subsection B, our choice of transformation [/ ( P, R) + XU is somewhat arbitrary, however it is 
sufficient to perform the exact diagonalization as shown in the present paragraph. Actually, since J7 (P, R) diagonalizes 
Hq at zeroth order in h, one can consider an-unknown- true diagonalization unitary operator Ui (P,R) reducing to 
?7(P,R) at zeroth order in h. That is Ui (P,R) = U {P,R) + n {...). As a consequence Ui (P,R) and C/ (P, R) +Xf/ 
are equal at the zeroth order in h and the difference 

6XU = Ui(P,Ii)-U(P,IL)-XU (30) 

is of order h. Moreover C/i (P, R) and ?7 (P, R) + XU being both unitary, SX is easily seen to be antihermitian. As 
a direct consequence, one can check that the difference between the exact diagonalization and our approximate one 

[{U (P, R) + XU) Ho (P, R) {U+ (P, R) + U+X+)] - [Ui (P, R) Ho ( P, R) U+ (P, R)] (31) 

is equal to 



[6X,e{P,R)] 



(32) 
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Given that e (P, R) is diagonal, this last term is always non diagonal. As a consequence, if we project on the diagonal 
the difference between our transformations, one gets 0. 

P [([/ (P, R) + XU) Ho (P, R) ([/+ (P, R) + U+X+)] - P [Ui (P, R) Ho (P, R) U+ (P, R)] = (33) 

Here we have denoted P[...] the projection on the diagonal. 

Now, given that Ui (P, R) Ho (P, R) (P, R) is truly diagonal, one has 

P [Ui (P, R) Ho (P, R) U+ (P, R)] - Ui (P, R) Ho (P, R) U+ (P, R) 



(34) 

so that ultimately 

P [{U (P, R) + XU) Ho (P, R) (C/+ (P, R) + U+X+)] ^ Ui (P, R) Ho (P, R) U+ (P, R) (35) 

We can therefore conclude, that the projection of our quasi-diagonalized Hamiltonian, by eliminating thus the non 
diagonal parts, is in fact the genuine diagonalized Hamiltonian at order h. 

F. The Diagonal Hamiltonian 

From the previous discussion we understand that the genuine semiclassical diagonal Hamiltonian Hd is simply 
given by the projection on the diagonal of the Hamiltonian Eql^D 

Hd = P[iUiP,K)+XU)HoiP,R){U+ iP,K) + U+X+)] 

= £ (P, R) + i [An,VH,e (P, R) + V^^e (P, R) + l IM^P^e (P, R) + Vp,£ (P, R) Ap,] 



-P 



-—[eiP,Il),Ap^AR, 



2h 



[e{P,-R),AR,] Ap, 



(36) 



where we introduced the notation A ^ P [A] . This Hamiltonian can be rewritten 

Hd = eiP+Ap,Tl+Ap) + ^P[[e{P,-R),AR,]Ap^[eiP,R),Ap,] Ap,] 

- e{p,r) + ^P[[eip,r),Ap,]Ap,-[eip,r),Ap,]Ap,] (37) 

where we have defined the projected dynamical operators 

r = P[(C/(P,R))RC/+(P,R)] =R+^R 

p = P[C/(P,R)P[/+(P,R)] = P + ylp (38) 
The non-canonical dynamical variables (p, r) have corrections of order h through the presence of the Berry connections. 



G. The equations of motion 

Given the Hamiltonian derived in the previous subsection, the equations of motion can now be easily derived. As 
usual [sl [lH the dynamics has to be considered, not for the usual position R and momentum P, but rather for the 
projected variables r and p. These new dynamical operators which naturally appear in our diagonalization process 
at the h order have components which do not commute any more. Actually 



[n,r,] = zQ^J = {Vp^Ap^ - Wp^Ap^) + [Ap^,Ap^] 

[p„Pj] = iQPP = -ih {Vp^Ap^ - Vp^Ap^) + [Ap^,Ap^] 

[p^,rj] = -ihd^j + iQYj = -i^^^3 - (Vfl.^fl, + Vp^^p,) + [Ap^,Ap^] 



(39) 



the Qij being the so called non-Abelian Berry curvatures. They are of order h as lowest order corrections to the 
commutators, but will actually induce semiclassical corrections of order h to the equations of motion. Indeed the one 
band-Hamiltonian Eq. [37] yields directly to general equations of motion for r, p : 



r = l[r,e(p,r)] + l 
p = l[p,e(p,r)] + l 



r, 2^^'[[£ (p,r),^pj^pi - [e (p, r) , ^pi] ^p,] 



P'^P [[£ (P, r) , Ap,] Api - [e (p+, r+) , Api] Ap,] 



(40) 
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where the commutators can be computed through the previous commutation rules Eql39l The last term in each 
equation represents a contribution of "magnetization" type (see the following applications) and has the advantage to 
present this general form whatever the system initially considered. However, to put some flesh on these equations, we 
now turn to several examples covered by our formalism. 

III. APPLICATION 1: THE DIRAC ELECTRON. 

To apply our previous formalism, we will consider two cases of Dirac Hamiltonians : The electromagnetic field and 
the static symmetrical gravitational field. These two cases have already been treated by different methods (111 and 
[13, [3), but in the second case (gravitational field) references to Berry phases was made for the first time in |l6|. 

A. The Dirac electron in an electromagnetic field 

The diagonalization of the Dirac Hamiltonian in the presence of an electromagnetic field is a difficult problem 
which was solved only approximately in the non-relativistic limit in an expansion. Another approach consists in 
diagonalizing the Hamiltonian at the semi-classical order as was done in [3| using an approximate Foldy-Wouthuysen 
transformation (19] . From the semi-classical Hamiltonian the equations of motion were derived showing a topological 
spin-transport effect due to the presence of the Berry phases. Here we propose to apply our general formalism for the 
semi-classical diagonalization of the Dirac Hamiltonian to show the effectiveness of our general method. It is worth 
noticing that the method developed above has now to be adapted as we will transform the Dirac Hamiltonian into a 
(2 X 2) block-diagonal matrix (due to the spin degree of freedom). However, it is easy to check that, since the non 
diagonal components for the energy blocks are of order 7i, these corrections do not impair our general formulas for 
the Berry phases contributions. 

We thus start with the following Dirac Hamiltonian: 

ffo (P, R) = a. (P-A(R))-h /3m + y(R) (41) 
where the matrix a and f3 are the usual (4 x 4) Dirac matrices: 




(42) 



with (Ti the usual (2 x 2) Pauli matrices (i = 1, 2, 3) and 



/ 1 
10 
0-10 

Vo -1 



(43) 



Replacing R by the parameter r, A(r) just shifts the momentum, so that we can diagonalize Hq (P,r) through the 
well known Foldy-Wouthuysen transformation 

t/(P,?) = ^±^^^±a=£=A(£)) (44) 
^2E{E + m) 



wheieE — y (F — A(r))'^ + m^. In this context, introducing the dependence in R we define the (non projected) Berry 
connections at first order in ?i as a (4 x 4) matrix 

^ la. (P - A(R)) P,/3 + iPEjE + m)a, - g (S x (P - A(R))), 
^^^-'^'US/p^U =n 2E^{E + m) ^^^^ 

and 

Ap, = -ihUVp, U+ = Vfl, Ak (R)Ar, (46) 
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where E will now denote the symmetrized form oi E = y (P— A(R))^ + m? (here at zeroth order in fi we can consider 
that P and R commute). The spin matrices (i — 1, 2, 3) are given by 

/ 

' 


V 

The general method developed in the previous section allows us to write the diagonal Hamiltonian iJ/j as a matrix 
generalization of Eg [571 

Hd = P [UH„U+] = Pe (p, r) + ^P[[e (p, r) , Ar,] Ap^ - [e (p, r) , Ap.] Ap,] + ^(r) (47) 



with £ (p, r) = Y (p — A(r))^ + vn?. In Eq. |47]the operators (p, r) are the physical dynamical variables satisfying the 
non-canonical commutations relations Eq. [39] Using the expressions for the Berry connections, Hp) can be rewritten 
as: 



[e (p,r),C/Vp.C/+] £*^''=t/Vp.C/+ 



-4^ + V{v) (48) 



Moreover, a straightforward computation shows that one can write P\_—^ [s (p, r) , UV piU'^] .£*-''^C/Vp^ JJ+j 
f3 ^^^ — at the first order in Ti, where we have introduced the intrinsic angular momentum of semiclassical 

particles L = Px^r, with = P[Ap] — fi ^'^2E(E+m)^ ^^'^ projection of the Berry connection on the diagonal. As 
a consequence, the Hamiltonian to be considered is given by 

, , nS.B L.B , , 

HD=pe{v.v)+(3^^-(3— + V{v) (49) 

which is the Hamiltonian deduced in [i*] from a different approach and which leads of course to the dynamics described 
in that paper. 

B. The electron in a static gravitational field 

The behavior of Dirac particles in static gravitational field is an important issue, at the crossroad of particle physics 
and cosmology. Different approaches for the diagonalization of the Hamiltonian leads to contradictory results in 
particular with regard to the existence of a dipole spin-gravity coupling • It is not our goal to discuss this 

specific point but we study the semiclassical diagonalization of the Hamiltonian to get the velocity and momentum 
evolution. We can in particular compare our results with the article [18] who uses a Foldy-Wouthuysen transformation 
instead of a semiclassical approximation. 

The interaction of a Dirac particle with a symmetric static gravitational field (500 = ^(R-);5io — ^T9ij — SijFCR)) 
is described by the Hamiltonian [18,] 

Hq = ^ (a.PF(R) + F(R)a.P) + /3mV{r) (50) 
The Foldy Whouthuysen transformation when R is replaced by a parameter r is given by 

y/2E{E + mV{r)) 



with E (P,r) — ■\/F2(r)P2 +m?V'^{r). This is quite the same as the free particle transformation. As a consequence 
introducing again the R dependence yields the non projected Berry connections for the position and the momentum 
operators : 

, ....^ ^ tF^R)a.PP,P + tpF{K)s{e + mV{R))a, - sF^jK) (S x P)^ 
Ar^ = inUVp^U — n 



2£2 (e + mV(R)) 

Ap. = -^nUS^P-U- = _,^r.(Vp.0)/3F2(R)a.P 

2e'^ 
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with (j) = ^, S the spin of the electron and the energy e is now given by e(P,R) = 

The expressions for Ar^ and Ap^ allow us ultimately to define the semiclassical transformation : J7 (P, R) + 
[Api , Ab^,] L/ (P, R) and to compute the diagonal Hamiltonian 

Hd = P [UHoU+] = (3e (p, r) - p^^mhV^{r). (p x E) (53) 

where the dynamical variables which are deduced from the projections vIr = -P[.Ar] = ^ 2£(e+mV') ^^'^ ^ -f[-^p] = 
are 

. . F2(R)S X P 

r = k+Ar = n n- ^ ' 



2e(£ + TOF(R)) 

p = P+v4p = P (54) 



The commutators between these variables are thus 



with the Berry curvatures 



= (55) 



2£3(p^r) \ e (p, r) + mV^(r) 

e^f = (56) 

One can check, after developing r as a function of R and the Berry connection, that the Hamiltonian Eq. [53] coincides 
with the one given in at order h. This confirms also the validity of the Foldy Wouthuysen approach asserted in 
in opposition with the transformation proposed in [l3|- However our approach is more general since it does not 
require an expansion in V and F as done in [isj . Of course, we retrieve the result of if we expand expression 1561 
at the leading order in F and V . Note also that when m = one recovers the Hamiltonian for the Neutrino or the 
photon proposed in [3j[l6|. 

To conclude this paragraph, we can derive the equations of motion with the help the noncanonical commutators 
between the coordinates and the spin 



(P,r) (£(p,r) + my(r)) 



eP^ = b.,I],] = (57) 

We then can deduce the semiclassical equations of motion for the electron in a symmetric static gravitational field 
(by projecting on the positive energy subspace) 

(AVI- \ 1 7 

1-- Vpe--pxe- + -V.£.0- 



VrE (58) 



where we have defined the vectors ©'"'" and 0''^ through 0[J = e^J'^Q'^^ and OfJ = e^^^Q]^ ■ 

The velocity equation in Eg 1581 contains in particular an anomalous velocity term pxQ^^ /h of order h which causes 
an additional displacement of the electrons orthogonally to the momentum p. The dynamics of the system must be 
completed by the spin dynamics which is 

ht = ^^M^sx (Vy(r)xP) - -Ex (V^^(r)xP) (59) 
e(e + mV[r)) e 

Although the position dynamics differs from the one of obtained in [l^, due to our choice for the physical position 
operator (r instead of R), one can show that the equations for p and hT, reduce to [Uj in the case of weak fields. 
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IV. APPLICATION 2 : THE ELECTRON IN A PERIODIC POTENTIAL 



This application has already been independently studied in [15| . and is easily recovered by the present general setup. 
The purpose is to find the semiclassical Hamiltonian for an electron in a periodic potential facing an electromagnetic 
field. This topic was also already dealt with in ^] in the context of wave packets dynamics. We will show that the 
semiclassical equations of motions which are very essential in the solid state physics context must be corrected by 
Berry phases terms. To apply our formalism, consider an electron in an crystal lattice perturbated by the presence of 
an external electromagnetic field. As is usual, we express the total magnetic field as the sum of a constant field B and 
small nonuniform part (5B(R). The Schrbdinger equation reads {Hq — e0(R)) 'I'(R) = i?5'(R) with Hq the magnetic 
contribution (cf) being the electric potential) which reads 



Ha=[— + eA(R) + e5A(R) 
' 2m 



P = -i?iV 



(60) 



where A(R) and (5A(R) are the vectors potential of the homogeneous and inhomogeneous magnetic field, respectively, 
and 1^(R) the periodic potential. The large constant part B is chosen such that the magnetic flux through a unit 
cell is a rational fraction of the flux quantum h/e. The advantage of such a decomposition is that for (5A(R) = the 
magnetic translation operators T(b) = exp(iK.b) defined in Appendix 2, with K the generator of translation, are 
commuting quantities allowing to exactly diagonalize the Hamiltonian and to treat 5A(R) as a small perturbation. 
The state space of the Bloch electron is spanned by the basis vectors of plane waves \ti, k) = |k) (g) \n) with n 
corresponding to a band index. The state |n) can be seen as a canonical base vector \n) = (0...010...0...) (with 1 at 
the nth position) such that U'^ (k) \n) — |u„ (k)) with |m„ (k)) the periodic part of the magnetic Bloch waves. In 
this representation K |n,k) — k|rt, k) and consequently the position operator is R —id/dii, implying the canonical 
commutation relations [Ri,Kj] —iSij. 

We first perform the diagonalization of the Hamiltonian in Ea l60l for SA = by diagonalizing simultaneously Hq 
and the magnetic translation operators T. The diagonalization is performed as follows: start with an arbitrary basis 
of eigenvectors of T. As explained in Appendix 2, in this basis Hq can be seen as a square matrix with operators 
entries. -ffo is diagonalized through a unitary matrix U (K) which should depend only on K (since U should leave 
K invariant, i.e., = K) and whose precise expression is not necessary for the derivation of the equations of 

motion, such that UHU^ ~ ^(K) — e(/<([/RC/+), where f (K) is the diagonal energy matrix made of elements £„(K) 
with n the band index (i.e. the diagonal representation of Hq). 

Now, to add a perturbation 5A(R) as in ([iBl)) that breaks the translational symmetry, we have to replace K in all 
expressions by 



K = K 



M(R) 



(61) 



and as the flux (5B on a plaquette is not a rational multiple of the flux quantum, we cannot diagonalize simultaneously 
its components Ki since they do not commute anymore. Actually 



(62) 



As a consequence of this non-commutativity, we just aim at quasi-diagonalizing our Hamiltonian at the semiclassical 
order (with accuracy fi). To do that we replace C/(K) by U (^K^ , so that the non projected Berry connections are 
Ar^ = iUV ^ and Aki = ^ Ri^AkCIVjAR^. From these we can define the nth intraband position and momentum 

operators r„= H+An and kn ~ K — eAn0^n)x i5B(r„)/?i + 0{h) with An = Pn{UV^U~^) the projection of the 
Berry connection on the chosen nth Band [l5|. It can be readily seen that the matrix elements of An can be 
written An (k) = i (u„ (k)| Vk |u„ (k)) (see also ref. HHfor the derivation of the position operator in the diagonal 
representation) . What is totally new here is the transformation on the momentum operator fc„ which get also a Berry 
connection correction. 

Using our results of section 2, the full Hamiltonian Eq. [60] can thus be diagonalized through the transformation 
C/(K) + jj- [Ari,Api] U(K) plus a projection on the chosen n-th Band as it is usual in solid state physics (the so 
called one band approximation) 



U 



HU+ K 



= Pn 



£ (k) - ^ [£{K),UVk,U+] e^i^^-^UVK,U+ 



-UVk,U+ [£{K),U\7K,U+]e 



ijk 



MiK).SBirn)+0{h^) 



(63) 
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where the energy levels 5„(k„J are the same as £'„(K) with k„ replacing K. The magnetization 7W(K) — 



We mention that this magnetization (the orbital magnetic moment of Bloch electrons) , has been obtained previously 
in the context of electron wave packets dynamics Q . 

From the expression of the energy Eq. [53] we can deduce the equations of motion (with the band index n now 
omitted) 



Some recent applications of semi classical methods to several branches of Physics, such as spintronics or solid 
state physics have shown the relevance of Berry Phases contributions to the dynamics of a system. However, these 
progresses called for a rigorous Hamiltonian treatment that would allow for deriving naturally the role of the Berry 
phase. 

This paper has been devoted to derive a semi classical diagonalization method for a broad class of quantum systems, 
including the electron in a periodic potential and the Dirac Hamiltonian. Doing so, we have exhibited a general 
pattern for this class of systems implying the role of the Berry phases both for the position and the momentum. In 
such a context, the coordinates and momenta algebra are no longer commutative, and the dynamical equations for 
these variables directly include the influence of Berry phases through the parameters of noncommutativity (Berry 
curvatures) and through an abstract magnetization term. Applications of our formalism consider the Dirac electron 
in an electromagnetic field, or in a particular case of static gravitational field, as well as the electron in a periodic 
potential. Our results are promising and indicate that our method will probably apply to several other systems. 

Acknowledgment. The authors wish to thank Aileen Lotz for a critical reading of the manuscript. 



r = dE{k)/hdk-kx e(k) 

hk = -eE - er X SB{r) - MdSB/dr (64) 

where [r*,r^] = i8'-'(k) with 6*^(k) — d''A^{k) — d^A'^{k) the Berry curvature. 

As explained in 15] these equations are the same as the one derived in [o] from a completely different formalism. 



V. 



CONCLUSION 



A. Appendix 1 



To start our computation, we will use the formula given in the preliminary : 



-;iAsymVfl,Vpz [C/(P,R)ffo(P,R)(7+(P,R)] = -UAsymV r,\/ pi [U (. 
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Let us first remark that iJolPjR-) is already symmetrized in R and P. As a eonsequenee 
AsymV pi [Hq (P,R)] = 0. Actually, remember that the asymmetrization term was the sum of com- 
mutators obtained by pushing the momentum half on the left and half on the right. For the same reason 

AsymVR, W pi [U (P, R)] = Asym^R, Vp, [U+ (P, R)] = (66) 

Now, we introduce the transformed variables and the non projected Berry Phases at order h: 

r = (C/(P,R) + X)R(C/+(P,R)+X+) ~R+ [i;i[/(P,R)VpJ7+(P,R)] =R + ^fl 

p = U{P, R) PC/+(P, R) ~ P - [iUU (P, R) VflC/+ (P, R)] = P + Ap (67) 

Before going further, we can find some relations on the Berry Phases. Given that : 

C/(P,R)?7+(P,R) = 1 (68) 

at the zeroth order in h we have the following relations at the first order in h : 

Ar = inU{P.Il)WpU+(P,Ti) = -iVpU{P,Il)U+{P,Il) 

Ap = -ihU{P,Ii)\7RU+(P,R)^iVR,U{P,K)U+{P,Il) (69) 

Using these results as well as V pi = ^ [Ri,], Vr, = ^ [P\] and inserting the operators U {P,TL) and 
(P, R) when needed, we get : 

'-hAsymVR,Vpi [U (P, R) Ho (P, R) U+ (P, R)] 
= ^Api [/ (P, R) [Ri, Ho (P, R)] U+ (P, R) + ^Ar, U{P,R) [Pi , Ho (P, R)] U+ (P, R) 

-^U (P, R) [Ri, Ho (P, R)] U+ (P, R) Api + (P, R) [P', Ho (P, R)] U+ (P, R) Ar, 

-^Apie {P,R) Ar, + ^AR,e {P,R) Api 
= ^Api [n,s (P, R)] + ^Ar, [p\ e (P, R)] - ^ [n,s (P, R)] Api + ^ [p\ e (P, R)] Ar, 

-^Apie{P,R)AR, + ^AR,e{P,R)Api 
= I i-^Ri^R'S (P, R) + V^ze (P, R) Ar,] + i [ApiVp,e (P, R) + Vp,e (P, R) Api] 

-^^p' [-Ar, , e (P, R)] - ^ [Ar, , e (P, R)] Ap, + ^Ar, [Ar, , e (P, R)] + 

^ [Apt , e (P, R)] Ar, - ^Ap,e{P, R) Ar, + ^Ar,s (P, R) Api (70) 
Rearranging the commutators leads to : 

^hAsymWR,Wp, [U (P, R) Ho (P, R) U+ (P, R)] 
= I [^K, Vij,£ (P, R) + Vii,£ (P, R) Ar,] + ^ [Ap,Vp,e (P, R) + Vp,e (P, R) Ap,] 

+^[£(P,R),^flJ^pi -^^pi£(P,R)^ii, + ^^ii,£(P,R)^pi (71) 

rewrite ^^p,e(P,R).AK, + ^Ar,s {P ,R) Api as -^.Ap; [e (P, R) , .AijJ + j^AR,[e {P,R) , Api] + 
jf^[AR„ Api] e{P,R) 

so that : 

'-nAsymVp, V pi [U (P, R) Ho (P, R) U+ (P, R)] 

J [AR,WRie (P, R) + Vp!£ (P, R) Ar,] + ^ [Ap; Vp,£ (P, R) + Vp,e (P, R) Api] 



2 



■^[e(P,R),Ap,]AK, +^[e(P,R),AfiJAp, + ^[Ap,,Ap,]e(P,R) (72) 
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And we thus have 



C/(P,R)Fo(P,R) U+ (P,R) 
= £ (P, R) + ^ [AR,VRie{P, R) + VRie{P, R) Ar,] + ^ [^p; Vp,£ (P, R) + Vp,e (P, R) Api] 

-^[£(P,R),^pz]^K, +^[e(P,R),^H,]^p, +^[^j{,,^p,]e(P,R) (73) 

as claimed in the text. 

VI. APPENDIX 2. HOW SOLID STATES PHYSICS FITS IN OUR FRAMEWORK. 

In solid state Physics, we assume that the Hamiltonian is invariant through a discrete group of translations, for 
example a group of lattice translations, whose elements have the form 

r(b)=exp(-^/ A,(R + r)dr 1 exp f^P.b) (74) 



where b is an arbitrary lattice vector. The eigenvalues of this operator are degenerated and have the form 

exp (ik.b) (75) 

where k belongs to some reduced dual lattice (a fraction of the dual lattice, i.e. a plaquette in solid state physics). 
We aim at defining the generators K of these transformations as 

K.b = log(r(b)) (76) 

so that we can define 

Ki = db,-K.h (77) 

However, this logarithm cannot be defined uniquely, and we will build an explicit choice in the sequel. 

To do so, we work with the extended representation, so that the considered state space is defined, similarly to the 
Dirac case, by : L'^{R^)^E, where L'^{B?) is seen as the set of functions of the variable k, running on B?. S is a vector 
space of infinite size representing the bands. For each value of k, T(b) is diagonal in with eigenvalues exp (ik.b). 
We define Ki as acting diagonally as the multiplication by ki. In other words, we have defined the momentum through 
an extension : K = K ^ Ids- Consequently the position operator R is acting as «Vk- Therefore the state space of 
the Bloch electron is spanned by the basis vector of plane waves |n, k) = |k} ® \n) with n corresponding to a band 
index. The state \n) can be seen as a canonical base vector \n) = (0...010...0...) (with 1 at the nth position). 

Turning now to the diagonalization process for the Hamiltonian, this last one can be performed independently for 
each value of k, since the Hamiltonian commutes with the translations. We can thus see the Hamiltonian as a set of 
square matrices indexed by k, each of them acting on each copy of E. As a consequence the diagonalization matrix 
is a Block acting on E for each value of k. 

Note that this diagonalization matrix can of course be seen as an operator U (K, R) or, and this is the point of 
view we adopt here, as a matrix acting on each copy of E, that is, a matrix U(K), whose entries depend on K only. 
Actually, the dependence in R appears in the non diagonal elements, and we can discard them if we consider this 
"half matrix, half operator" version. This mixed representation has the advantage to do the connection with the 
Dirac Hamiltonian. 

In this set up, K being diagonal and proportional to the identity, it commutes with every matrix f/(K, R) 
preserving the Blocks. When considering the diagonalized Hamiltonian e(R, K), it can also be seen as a diagonal 
matrix (implicitly denoted e(K)) whose components are diagonal and denoted £„(K), the n th band energy. The 
commutator 

[R,£(K)] = V£(K) (78) 

is again a diagonal matrix whose entries are V£ri(K) (sketch of proof : £(K) is a series whose elements are products 
of powers of K and R. For each power of K, the fact that £(K) is diagonal implies that the dependence in R is 
a diagonal matrix. The gradient in K acting only on the power of K, this diagonality is preserved.) Turning now 
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to the perturbation SA^R), let us remark that if the operator (5A(R) preserves the bands, all operators e{'K)SA{IV), 
[e{'K),dA{IV)] are diagonal, and given our previous remarks, the same is true for Ve(K)SA{IV), [Ve(K),(5A(R)]. 

All this remarks that are obviously true for the Dirac case, appear to be useful in the solid state physics case 
(application 2), since it shows that the electron in an periodic potential fits in our framework. Actually, we can 
formally consider the Hamiltonian of such a problem as given by a matrix depending on the momenta and the 
coordinates. 
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